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Abstract 

This paper asks what classes of input signals are sufficient in order to completely identify 
the input /output behavior of generic bilinear systems. The main results are that step inputs 
are not sufficient, nor arc single pulses, but the family of all pulses (of a fixed amplitude 
but varying widths) do suffice for identification. 
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1 Introduction 



In this paper, we address the following question: what types of input signals are sufficient 
to completely identify the i/o behavior of a system? In other words, we look for classes IA of 
inputs with the property that, if a system a is stimulated with the inputs from the set hi and the 
corresponding time record of outputs is recorded, then, on the basis of the collected information 
on inputs and outputs it is possible — at least theoretically, with no regard to computational 
effort, and in the absence of noise — to obtain a system a which is equivalent to a (Figure 
By "equivalent" we mean that the estimated system a will be completely indistinguishable from 
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the true system a in its i/o behavior, even when presented with inputs that do not belong to the 
restricted class U used for the identification experiments. Whether a certain class IA of inputs 
is rich enough for identification is heavily dependent, of course, on prior assumptions about the 
system a. 

It is often very difficult to perform experiments in which arbitrary input profiles are used. 
Often, the only possible experiments are those in which steps, i.e. constant inputs, are applied. 
For example, in molecular biology, a step input corresponds to subjecting a cell culture to 
a fixed concentration of an extracelular ligand such as a drug or growth factor. Sometimes 
somewhat more complex inputs, such as pulses (keep the input constant at some level, then 
change it back to some default value) can be used, but this is already not easy to implement, 
much less more complicated test signals. This presents a theoretical challenge: how does one 
know if all possible "identifiable" information about the system can be obtained from such a 
restricted class of experiments? In the case of linear models of systems, this issue does not 
arise, because basically any single input, as long as it is nonzero, for example a single step or 
a single pulse, suffices for identification (or several inputs, one for each input channel, if the 
system has multiple inputs; for simplicity we restrict here to single-input systems). Note that 
we are thinking here of an experimental setup in which observations are collected over time. 
If, instead, only steady-state behavior was observed, and not transients, then one input is not 
enough, even for single-input linear systems. In that case, one has to use multiple inputs, such 
as steady-state measurements of responses to periodic signals at different frequencies. 

For nonlinear systems, it is thus an interesting question whether constant inputs or pulse 
inputs, or simple combinations of these input classes, suffice for identification, as they do for 
linear systems. In this paper, we show that for a large and interesting class of nonlinear systems, 
that of bilinear systems, constant inputs do not suffice, but pulses do. 

Bilinear systems constitute an appealing class of nonlinear systems [31 El EES El- While for 
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linear systems the evolution of the states is only allowed to depend on linear functions of the 
state variables and inputs, in bilinear systems one also allows a linear dependence on products 
between input and state variables. Bilinear systems can be easily described in linear-algebraic 
terms, and a theory, in many respects analogous to the linear theory, can be developed for the 
analysis of their input /output properties. On the other hand, bilinear systems are theoretically 
capable of approximating arbitrary input/output behaviors on finite time-intervals |211 IU1 ITH]. 
They have been used to model chemical processes, electrical networks, power plants, nuclear 
reactors, robotic manipulators, and many other systems in engineering, chemistry, biology, 
economics, and other fields |TI|. They can also be employed in order to model and analyze 
certain simple enzymatic signaling cascades, when substrates are not too close to saturation 
and thus Michaelis-Menten kinetics can be replaced by bilinear expressions jlO| I2f)j. 

Informally (see next section for precise definitions and statements), the main results that 
we prove are as follows. On the negative side, we show that step inputs are not enough for 
identifying bilinear systems, nor do single pulses suffice. On the positive side, we show that the 
family of all pulses (of a fixed amplitude but varying widths) do suffice for identification. To 
be precise, one must impose certain non-degeneracy conditions on the classes of systems being 
considered, and, for the negative result, one wants to avoid trivial counter-examples in order to 
say something interesting. Thus, all results are stated for generic classes of systems. 

Our techniques are based on realization theory. We make heavy use of ideas originally 
developed by Kalman for realizations by linear systems, and refined by Isidori and later Fliess 
for bilinear systems. 

The organization of this paper is as follows. Section |2] gives the basic definitions, and 
reviews the background from bilinear realization theory. We provide a self-contained discussion 
because, even though the results proved in that section are not new, it is hard to find references 
presented as needed here. The main results are stated in Sectional The proofs of the negative 
results are given in Section 0] and the proofs of the positive results are given in Sectional The 
latter are inspired by Juang's very nice paper (we thank David Angeli for bringing this 
paper to our attention). Although the bilinear identification algorithm given in ^2] involves 
some ambiguous steps, such as taking (non-unique) logarithms of matrices, we were able to 
adapt many of its basic ideas; we discuss in Section |3 connections to that work. Conclusions 
and some remarks are presented in Section [Tj 

2 Preliminaries 

Systems and I/O Equivalence 

Our results will be for bilinear systems, which are defined by affine vector fields on M. n and hence 
are described by matrices, but an abstract setup will allow us to discuss some preliminary facts 
in more generality. We consider single-input single-output initialized systems a, in the usual 
sense of control theory |19j : 

x(t) = f (x(t)) + «(t)/i(a:(t)), x(0) = x , y(t) = h(x(t)) (1) 

(we will drop the arguments "(<)" if clear from the context), where /o and f\ are smooth 
vector fields on a manifold M and h is an output function M — > R. Inputs can be taken to 
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be any Lebesgue-measurable, essentially bounded, functions u : [0, T u ] — > M, but there will be 
no difference in results if one restricts to, for instance, piecewise continuous inputs. We let Q, 
be the set of all inputs. In principle, solutions are unique but only defined on some maximal 
interval (which depends on the initial condition and the input), but for simplicity, and since 
it is the case anyway for bilinear systems, we assume that solutions are defined for all times 
t G [0, T u ] (forward completeness). We let tp(t,u), or just x(t) if the input is clear from the 
context, be the solution of Q at time t, and y(t) = h{ip{t, u)) the corresponding output. When 
more than one system is being considered, we use appropriate notations; for example, a system 
a will be described by equations x = fo(x) + ufi(x), x(0) = x~q, y = h(x). 

Given two systems a, a, and an input u, we say that a, a are indistinguishable under u if 
h(cp(t,u)) = h((p(t,u)) for all t G [0, T u ]. If h(np(t,u)) ^ h((p(t,u)) for some t, we say that u 
distinguishes among the two systems a, a; in other words, the "input /output experiment" con- 
sisting of perturbing the system with this input u results in a different time-varying observation 
for a than for a. 

Given a subset IA C of inputs, we say that the two systems a, a are input/output (i/o) 

equivalent with respect to all inputs in IA if no input in IA is able to distinguish between the two 

systems a, a, i.e., provided that h(<p(t,u)) = h((p(t,u)) for each u £W and for each t G [0, T u ]. 

We write a = a in that case. 
u 

In the special case in which U = f2, we write simply a = a and simply say that the two 
systems are i/o equivalent. That is to say, they cannot be distinguished in any way whatsoever 
based on their "back box" input/output behavior. 

Let £ be a class of systems. A subset IA C of inputs is said to be sufficient for identifying 
systems in the class £ if, for any pair of systems a, a in £, 

a = a a = o . 
u 

In other words: whenever a is not i/o equivalent to a, there must exist some input in the set 
IA which distinguishes among the two systems a, a. 

Linear systems (finite-dimensional, continuous-time) are those for which /q is linear, fi is 
constant, and h is linear, i.e. systems described by equations 

x = Ax + bu , x(0)=0, y = cx (2) 

with A G ]R nxn , b G ]R nxl , and c G IR lxn . We often refer interchangeably to a linear system 
or its corresponding triple of matrices (A, b, c). Linear systems can be identified by any single 
nonzero input on a nontrivial interval, such as a constant function (a step), or a pulse. 

4- Tuples and Bilinear Systems 

We consider two slightly different classes of bilinear systems. To define these classes, we first 
introduce J^-tuples as follows: 

(A, N, b, c) where A, N G K nxn , b G K nxl , c G E lxn (3) 

(the integer n is called the dimension of the 4-tuple). 
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We say that a system Q is a bilinear system of type I if /o is linear, f\ is affine, xq = 0, 
and /i is linear. In other words, the system equations are: 



x = (A + uN)x + bu, x(0) = 0, . . 

y = cx, 

where (A, N, b, c) is some 4-tuple as in ©. We use a notation such as "cr°" to refer to a system 
of type I. With some abuse of terminology, we also simply write a° = (A, N,b,c). Note that 
linear systems © constitute the subclass of bilinear systems of type I for which N = 0. 

We say that a system Q is a bilinear system of type II if /q and /2 are both linear and h 
is linear (but the initial state may be nonzero). In other words, the system equations are: 

x = {A + uN)x, x(0) = b, 
y = cx, 

where (A,N,b,c) is a 4-tuple as in (J3J). Once again, we do not differentiate between a system 
of type II and its associated 4-tuple when the meaning is clear from the context. 



Similarity 

We say that two 4-tuples (A, N, b, c) and (A, N, b, c) are similar if they have the same dimension 
n and there exists an invertible n x n matrix T such that the following equalities hold: 

A = TAT~ 1 , N = TNT~ 1 , b = Tb, c = cT~ 1 . (6) 

Note that, for linear systems (N = N = nX n), this reduces to the familiar equivalence relation 
in linear systems theory. 

We say that two bilinear systems a, a (both of type I or both of type II) are similar (or 
"internally equivalent" ) , and we write 

<7 ~ o 

if there is a change of variables x = Tz such that the equations of a get transformed into those 
of <7. For systems of type I, this means that 

(A + uN)x + bu = x = Tz = T[(A + uN)Tx + Tbu 

for all x and u, and also cx = cTz = cz; thus, a ~ a is the same as saying that the 4-tuples 
(A,N,b,c) and (A, N ,b,c) are similar. An analogous statement holds for systems of type II. 

An easy calculation shows that a ~ a => a = a, and a converse holds as well, under certain 
minimality assumptions, as discussed below. 



Checking I/O Equivalence 

For analytic systems, input/output equivalence can be verified by checking certain algebraic 
equalities, and there is no need to test all possible inputs, as we discuss next. 

For any smooth vector field / on W l , and any smooth function h : M. n — > R, the Lie 
derivative Lfh : R n — > R is defined as the function (Lfh)(x) = Vh(x) f(x), where V/i is the 
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gradient of h. (In differential-geometric terms, Lfh is simply the value of the vector field on 
h, when vector fields are viewed as derivations on spaces of smooth functions.) More generally, 
if f\, ■ ■ ■ , fk are vector fields, the iterated derivative Lf k ... Lf t h is defined recursively by the 
formula L fk {L flt _^ . . . L fl h). 

Suppose that a and a are two systems (^Q) for which the vector fields /o and f\ are analytic 
and the function h is also analytic. Then, a = a if and only if 

(Lf lk ■ ■ ■ L fn h) (xo) = (l^ . . . Lj l) (T ) (7) 

for all sequences (ii,...,ik) £ {0, l} fc and all k > 0. (When k = 0, (J7J) says that h(xo) = h(x~o).) 
This is true because the expressions in (J2J) are the coefficients of the Fliess generating series 
of the input/output behavior associated to the respective systems, and the i/o behavior is in 
one-to-one correspondence with the coefficients of the series, see |23) . Lemma 2.1. 

For bilinear systems of types I or II, i/o equivalence amounts to an equality of vectors. 
Indeed, take first systems of type I. In this case, fo(x) = Ax, fi(x) = Nx + b, and h(x) = ex. 
Therefore, one can see inductively that: 

(l L h) (x) = i if ik = ° : cAil ' ' ' AikX 

\ /**■•■ f h ){ x ) \ if z fc = 1 : cA h . . . A ik x + cA h . . . A ik _ t b 

where Aq = A and A\ = N. In particular, for x = xq = 0, we have that (Lf. . . . Lf i h)(x) = 
for all sequences with = 0, and [Lfc . . . Lft h)(x) = cA.- Ll . . . Ai k _ 1 b for all sequences with 
ik = I- 

Generally, given two 4-tuples (A, N, b, c) and (A, N, b,c), let us say that they are i/o equiv- 
alent if 

cA h ...A ik b = cA h ...A ik b (8) 

for all sequences of matrices Aj picked out of A and N, including the "empty" sequence (cb = 
cb). (It suffices to check sequences of length n + n, where n,n are the respective state-space 
dimensions; cf. U HH El 113 HE] •) 

Then the preceding discussion proves: 

Lemma 2.1 Two systems a and a of type I are i/o equivalent if and only if the corresponding 
4-tuples are i/o equivalent. □ 

For bilinear systems of type II, the same conclusion holds, in this case because 

( L h k --- L f n h)(b) = cA n ...A lk b. 

Lemma 2.2 Two systems a and a of type II are i/o equivalent if and only if the corresponding 
4-tuples are i/o equivalent. □ 
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Canonical Systems and Uniqueness 

A 4-tuple (A, N, b, c) as in © will be said to be canonical provided that the following two 
properties hold: 

1. There is no proper subspace of M. n that contains b and is invariant under x t— > Ax and 
x i — > Nx. 

2. There is no nonzero subspace of M. n that is contained in the nullspace of x i— > cx and is 
invariant under x t— > Ax and x i— > Nx. 

The first property can be equivalently expressed by saying that the set of vectors of the form 

A h ...A ik b, (9) 

ranging over all matrix products with Aj € {A,N} (including k = 0, i.e., b), or equivalently 
over all products of length k at most n — 1, must span all of W 1 . This property is often called 
"span-reachability" because, for bilinear systems, it amounts to the requirement that the set of 
states reachable from the origin span all of the state-space. Similarly, the second property can 
be equivalently expressed by the dual property that the span of the vectors 



(prime indicates transpose) be all of W 1 (once again, length < n — 1 suffices), and is an observ- 
ability property for bilinear systems. Canonical 4-tuples are also called "minimal" because [31 
l51 llll l51 ll7U18| they have minimal dimension among all other 4-tuples which are i/o equivalent 
in the sense of (jHJ); moreover, if a 4-tuple (A, N, b, c) is not canonical, then [31 El EH El El H!8] 
there is some 4-tuple (A, N, b,c) which is canonical and is so that @ holds. (We do not need 
in this paper the interpretations in terms of reachability and observability, nor the minimality 
result.) 

We will call a bilinear system a (of type I or II) canonical if the corresponding 4-tuple is 
canonical. 

A very special case is that of linear systems (J2J), i.e. systems of type I with N = 0). Such a 
system a = (A, b, c) is canonical if and only if it is reachable and observable in the usual sense 
of control theory JS]. The controllability matrix 7Z(A,b) and the observability matrix 0(A,c) 
are defined respectively by: 

K(A,b) = {b Ab ... A n ~ l b), 

0(A,c) = K{A',c'y = (c 1 A'c 1 ... {A'Y^c')' 

(prime indicates matrix transpose) . The system a is canonical iff both matrices have full rank 
n. 

Similarity and I/O Equivalence 

We already remarked that a ~ a => a = a for any two bilinear systems (both of the same 
type). Conversely, if both systems a and a are canonical, a = a => a ~ a. Thus: 

if a and a are canonical, a = a a ~ a . (10) 
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This is a standard fact about bilinear systems 03 El QZI E] (strictly speaking, these 
references deal with discrete-time systems such as x(t + 1) = (A + u(t)N)x(t), but the algebraic 
statement about 4-tuples is the same as in the continuous-time case). The proof is, in fact, 
completely analogous to the proof for linear systems F° r completeness, we provide a proof 
here: 

Lemma 2.3 Suppose that the two 4-tuples (A, N,b,c) and (A, N ,b,c) are canonical and i/o 
equivalent. Then they are similar. Moreover, the similarity transformation T in © is unique. 

Proof. Pick any x € M. n . By the span-reachability property @ for (A,N,b,c), there are real 
numbers X a , where a denotes sequences (ii,...,ik) of length at most n — \ (including the 
"empty" sequence) such that x = Yl a ^aA a b, where we denote A a = A^ . . . Ai k for a = 
(ii, . . . , ifc). Now define Tx := X a A a b. 

There are many possible representations of a vector x as a linear combination of the spanning 
set in (J§J) for (A, N, b,c), so to see that T is well-defined as a mapping we need to verify that if 
Y^ a X a A a b = Y2 a ^'a^ab then ^ a \ a A a b = ^ Q \' a A a b. By linearity, it is enough to show that 
J2 a X a A a b = => ^2 a \ a A a b = 0. Suppose that ^ Q X Q A a b = 0. Then also cAp ^ a X a A a b = 
for any other index /3, or equivalently X a cAp a b = 0, where (3a is the concatenation of the 
sequences j3 and a. Now, i/o equivalence of the two given 4-tuples implies that cAp a b = cAp a b 
for all indices, and so also cApx = ^ Q X a cAp a b = 0. This holds for any index f3, so, using the 
observability of (A, N, b, c), we conclude that x = 0, as desired. The mapping T is obviously 
linear (by definition), and it is onto because of the reachability of (A, N, b, c), which means that 
every x £ W 1 can be written as J2 a X a A a b for some A Q 's. To prove that T is one-to-one, we 
simply reverse the argument used to prove that T was well-defined. Uniqueness follows by the 
same argument. I 

By picking among all the possible linear combinations the one whose coefficients have min- 
imal Euclidean norm, one obtains an explicit expression for T: 

T = TZTZ*, 

where # denotes matrix pseudoinversion, 1Z is a matrix listing the products in (jHJ) of length 
< n — 1, and 1Z lists the vectors in the same order for the second 4-tuple. For linear systems 
a = (A, b, c), the equivalence becomes 

T = K{A,b)K(A,b)- 1 , (11) 

where R is the usual reachability matrix ( |19j. Theorem 27). 

Generic Sets of Systems 

We will make statements about "generic" classes of systems, so we must define this term 
carefully. Genericity can be defined in many ways, for example in probabilistic terms (a set is 
generic if it has "probability one" ) or, as usual in mathematics, in terms of open dense sets. In 
order to provide the strongest possible results, we combine both definitions and say here that 
a subset S of an Euclidean space M. e is generic provided that: 
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• the set S has full measure, that is, the complement S c has Lebesgue measure zero, and 

• the set S is open (and dense) in Mr. 

When dealing with sets of 4-tuples ©, we view such sets as subsets of Mr with t = 2n 2 + 2n. 

When talking about genericity of classes of systems of type I or II, we mean genericity of 
the sets of associated 4-tuples. Specifically, if we let be the class of n-dimensional bilinear 
systems of type I, then we think of as M? n +2n , and similarly for the class Sj/ of n-dimensional 
bilinear systems of type II. 

3 Statements of Main Results 

For any a€l and any r > 0, let Z T ,a denote the class of all functions of the form 

a for < t < t, 



, . (a for < t < 



where (3 is a constant. Let u TiQ , denote the particular pulse function in T T>a for which (3 = 0, 
that is, 

, s (a for < t < r, , . 

U ^ (t) = l0 fort>r. (12) 

Note that in the special case when r = 0, X TOi becomes the class of constant functions. (There 
is a small ambiguity in that we have not specified the domain of the inputs. We can view these 
inputs as defined on some interval [0,T] with T > r; any such T will give the same results.) 

Let us now state the negative main results of this paper. 

Theorem 1 For any r > and any a G M, there is a generic subset S of such that, for 
every system a° € 5, there is some a° G S such that 

1. a° and a° are i/o equivalent under the pulse function u Ta (o~° = a°), but 

{u r ,„} 

2. a° and a° are not i/o equivalent. 

Theorem 2 There is a generic subset Q of S^ 1 such that, for every system a G Q, there is 
some a G Q such that 

1. a and a are i/o equivalent under all the pulses in the set I T>a (a° = o~°), but 

2. a and a are not i/o equivalent. 

By setting r = for the collection X T)0! , one obtains the following as a consequence of 
Theorem |2 

Corollary 3.1 There is a generic subset Q of S^ 1 such that, for every system a G G, there is 
some a G Q such that 
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1. a and a are i/o equivalent under every constant input, but 

2. a and a are not i/o equivalent. □ 

The first part of Corollary 13.11 mav be restated as follows: for every (A,N,b,c) £ Q, there 
exists some (A, N, b,c) € Q such that 



ce 



(A+f3N)t b = - e (A+/3iV)t^ 



for all (5 £ R and all i > 0. This implies that: 

c f ^ A+(3N ^-^f3bds = c f ; e &+P")* fids. 
Jo Jo 

for all (5 £ M. and all t > 0. Hence, the result in Corollary 13 . II also applies to systems as in 

Corollary 3.2 There is a generic subset Q of such that, for every a £ Q, there is some 
o £ ^ such that 

1. <r and a are i/o equivalent under every constant input, but 

2. a and a are not i/o equivalent. □ 

Next, we state our positive results for systems of both types. For any a £ K, let V a denote 
the set of pulses of magnitude a: 

V a := {u TiQ | r > 0}. 

Theorem 3 For each a / 0, there is a generic subset M of such that, for every pair of 
systems a°, £ M., 

a° = a° a° = a° . 

v a 

Theorem 4 For each a / 0, there is a generic subset Ai of S*J such that, for every pair of 
systems o~°, a 2 £ M., 



4 Proofs of Negative Results 
4.1 Some Preliminaries 

The following construction is key to the proofs of the negative results. The following observation 
was apparently first made in [2j (see problem 1 in page 110, and problem 12 in page 105). 

Lemma 4.1 For each canonical triple a = (A,b,c), there is a unique matrix T = T{a) such 
that 

AT = TA' , b = Tc , cT = b' . (13) 
Moreover, the matrix T(a) is given by T(o~) = 7Z(A,b) [(0(A, c))'] _1 . 

10 



Proof. Observe that for each canonical triple (A, b, c), the triple (A, c' , b') is also canonical, and 
the two triples are i/o equivalent since cA k b = (cA k b)' = b'{A') k d for all nonnegative integers 
k. Thus there is a (unique) similarity between (A,b,c) and (A,c',b'), an invertible matrix T 
such that: 

AT = TA' , b = Tc' , cT = b' . 

The formula for T is given in which, since in this case A = A' and b = c' and 0(^4, c)' = 
lZ(A',d), reduces to that shown. I 

Although not needed, it is worth remarking that T is symmetric. This can be proved as 
follows: transposing the relations in (|13|) . one has that also AT' = T'A', b = T'c' , and cT' = b' . 
Since the similarity T is unique, T = T' . For each nonzero nxn matrix S, consider the following 
set: 

B(S) := {N € R rtxn | NS = SN'}. 

Note that B(S) is a proper linear subspace of M nxn , because, in particular, when Eij is the 
matrix having a 1 in its (i,j)th position and zero elsewhere, EijS = SE[ - implies that s^j = 
for all k ^ i. 

We now define a set Qq that will play a major role in the constructions. It is defined as the 
set consisting of those 4-tuples 

(A,N,b,c) € R nxn x M nxn x R nxl x M lxn 

such that: 

1. (A, b, c) is canonical, 

2. N &B(T(A,b,c)). 

Note that, since the triple (A, b, c) is already canonical, every element of Qq is canonical as 
a 4-tuple. 

For the next result, it is more elegant not to use inverses. Given a triple a = (A,b,c), let 
O denote the cofactor matrix of O = 0(A,c). Note that if a is observable, then O" 1 = j^O, 
where A(tr) = detO. 

For any triple a = (A,b,c) (not necessarily minimal), we define T(a) := 7Z(A,b) (0(A,c))'. 
This is a polynomial expression on the entries of A, b, and c. If (A, c) is observable, T(a) = 

Observe that B(S) = B{5S) for any scalar 5. Thus, B(T(A,b,c)) = B{f(A,b,c)). 

Lemma 4.2 The complement Qq of Q$ is a proper algebraic subset of M. nxn x M nxn x M nxl x 
M lxn . 

Proof. The complement of Qq of Qq is the union of the solution sets of the following equations 
respectively: 

det K(A, b) = 0, 
det 0{A,c) = 0, 
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and the n 2 scalar equations given by 

NT(a) = f(a)N'. (14) 

Hence, Gq is an algebraic set. Each subset is proper (for the last one, pick an arbitrary canonical 
(^4, 6, c) and refer to the above remark that B(S) is always proper), and hence of dimension less 
than n 2 + 2n, so the union is also proper. I 

Lemma 4.3 For each (A, N, b, c) G Go, consider M := TN'T" 1 , where T = T(A, b, c). Then, 

1. M^N, 

2. (A,M,b,c) G Go, and 

3. for each 7GI and each nonnegative integer k: 

c(A + jN) k b = c^A + ^Mfb. (15) 

Proof. Let (A, N, b, c) G Go- The fact that M^N follows from the fact that N B(T(A, b, c)). 
To see that (A, M, b, c) E G , note that MT = TN'T^T = TN' ± TM' since N ^ M and T is 
invertible. The equality (HSJ) follows by the equalities AT = TA', b = Tc 1 , cT = b', TN' = MT, 
and the following: 

c{A + -fN) k b = {c(A + ~fN) k b)' = b'(A' + -fN') k c' 

= cT (T~ 1 AT + -fT~ 1 MT) k T~ x b = c(A + -yM) k b 

for all k > and all 7 G R. I 

Corollary 4.4 For each (A,N,b,c) G Go, and M = TN'T" 1 , the 4-tuples (A, N, b, c) and 
(A,M,b,c) and are not i/o equivalent to each other. 

Proof. Suppose that these two 4-tuples would be i/o equivalent. By Lemma [2 .'61 they are similar. 
Let T provide a similarity as in (JSJ). In particular, T provides a similarity between the canonical 
triple (A, b, c) and itself. Since there is a unique such similarity, and the identity / is one, it 
follows that T = I. Thus N = TMT' 1 = M, contradicting the fact that N ^ M. I 

4.2 Proof of Theorem [T] 

Let C be the subset consisting of all those 4-tuples 

(Q,N,b ,c) G rn xr xn xr xl xR lxn 

which satisfy the following conditions: 
(a) (Q, 6q , c) is canonical, 
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(b) (Q — N, bo, c) is canonical, 

(c) — I is invertible, and 

(d) N<£B(T(Q,bo,c)). 

Letting A(F) denote the collection of eigenvalues of a matrix F, the Spectral Mapping 
Theorem implies that e A ^) = {e A : A 6 A(Q)}. Thus, assumption (c), which says that 1 is not 
an eigenvalue of e®, implies, in particular, that Q is invertible. 

Lemma 4.5 The complement C c of C is a countable union C c = (Jfc^=o^ fc ' wnere each E\~ is a 
proper algebraic subset of l ttXn x R ttXn x l nxl x R lxn . 

Proof. First note that C Q Go, and C c is the union of Gq and of the solution sets of of the 
following equations: 

detK(Q - N,b ) = 0, (16) 
det 0(Q - N, c) = 0, (17) 
det(e Q -J)=0. (18) 

Clearly, the solutions sets A\ and A2 of Equations (JT5|) and (|T7|) respectively are proper algebraic 
sets. By the Spectral Mapping Theorem, (|18|) holds if and only if — 4/c 2 7r 2 is an eigenvalue of Q 2 
for some integer k, and hence, the solution set A3 of (|TH|) is the countable union of the solution 
sets {^3fc} of the equations 

det(Q 2 + 4fcVl) = 0, k = 0,1,2, ... . 

Hence, C c is the countable union of Gq, Ai, A2, and {-4.3fc}fc>o- ■ 

Let X = R nxn x R nxn x R nxl x M lxn , and consider the analytic map ip : X -» X defined 

by 

tl> : (Q, N, b ,c) ^(Q-N, N, (p(Q))*b , c) , 

where p(Q) = Jq e s ®ds, and p(Q)* denotes the adjoint matrix of p(Q). Note that if — I is 
invertible, then Q is invertible, and the matrix p(Q) = Q~ 1 {e ( ^ — I) is also invertible. Hence, 
when restricted to the open set 

%i := {(Qi N, bo, c) : — I is invertible} , 

ip is given by 

(Q, N, bo,c)^(Q- N, N, det(p(Q)) (p{Q))-\, c) . (19) 
Let V'o denote the restriction of tp to X\, and consider the open set 

3>i := {(A, JV,6,c) : (e A+JV - J) is invertible} . 
Then, ipo is a (smooth) diffeomorphism from Ai to y±. Its inverse is : y% — ► A4 given by 

^0 1 : (A, iV, 6, c)~U + iV, JV, det e s(A+7V) ds^ * e s(A+7V) ds^ b, . (20) 
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Since y = ip{Xi) C ip(X), it follows that the complement (ip(X)) c of is a subset of 

{(A, iV, b, c) : (e^ +Ar — /) is singular}. Hence, (ip(X)) c is contained in the countable union of 
the solution sets jF fc of the equations 

det({A + N) 2 + 4k 2 ir 2 I) = 0, k = 0, 1, 2, . . . . 

Let P = ^(C), and write C c = UfcLo ^fc> wnere the sets £}~ are as in Lemma l4~51 The next lemma 
then follows from the fact that V c C [-0(^)] c (J -^(C c ), and Lemma 14.51 



Lemma 4.6 P c C (^U^Lo V'o(^fe)J U(^UfcLo-^J' w here ipo is the diffeomorphism from X\ to 
3^1 defined hv ljT3jl . □ 

Corollary 4.7 is generic. □ 

Proof. Since every proper algebraic set has measure zero, the set IJfc>i ^"k has measure zero. 
Furthermore, since the image of a measure zero set under a differentiable map has measure 
zero (see e.g. pQ, Lemma 2.6), ip(£k) has measure zero for each k. This implies that T>° has 
measure zero. Therefore, T> is of full measure, and as a consequence, T> is dense. Finally, T> 
is open because it is the image of C under the diffeomorphism -00 and 34 is an open subset of 

R nxn x p xn x R nxl x R lxn , 

Let u = u TjQ , with r = 1, a = 1. 

Lemma 4.8 Consider systems as in Q). For every cr° £ T>, there exists o~° £ V such that the 
following holds: 

1. a° and cj° are i/o equivalent under the pulse function u, but 

2. (j° and a° are not i/o equivalent. 

Proof. Let cj° = (A, N,b,c) £V = ip(C). Then there exists (Q, N, b , c) £ C C Go such that 

A = Q-N, b=[det(p(Q))](p(Q)r 1 b . 
Let b\ = p(Q)b. Then b\ = det(/9((5))6o, and hence, 

n{QM) = det(p(Q))K(Q,b ) 
TZ(Q-N,h) = det(p(Q))K(Q - N,b Q )- 
Since det(p(Q)) 7^ 0, both (Q,&i) and (Q — N,bi) are reachable. Moreover, 

T(QM,c) = T(Q,det(p(Q))b ,c) 

= W,det( /0 (Q))6 )[O(g,c)r 1 
= det(p(Q))n(Q,b )[O(Q,c)']- 1 
= det(p(Q))r(Q,6 , C ) 5 

which implies that B(T(Q,bi,c)) = B(T(Q,bo,c)). Therefore, (Q, N,bi,c) G C. In particular, 
(Q,N,b,c) e g - 
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Applying Lemma f-Pl to (Q, N, b\, c), one sees that with M = TN'T 1 (where T = T(Q, bi,c)), 
it holds that M ^ N, and 

ce *(Q+7A0 6l = ce t(Q+ 7 ^) 6l v 7 , V< > 0. 

In particular, for 7 = —1, 

ceW- N \ = ce^- M \ Vt>0. 

Let a° = {{A + N — M),M, b, c). Consider the two systems a and a°: 

x = (A + uN)x + bu, x(0) = 0, y = cx, 

z = [(A + N - M) + uM]z + 6n, z(0) = 0, y = az. 

Since u = 1 for t € [0, 1], the two systems reduce in that interval to: 

x = (A + N)x + bu, x(0) = 0, y = cx, 0<t<l 
z = (A + N)z + bu, z(0) = 0, y = cz, < t < 1. 

It follows that x(t) = z(t) on [0, 1], and hence outputs coincide for t G [0, 1]. In particular, at 
time t = 1 both systems are in state 

f 1 eQ-'W+Wbds = f 1 e^-'Wbds = p(Q)b = h . 
Jo Jo 

Now, for t > 1, using that A = Q - N and A + N - M = Q - M, we have that 

x = (Q-N)x, x(l)=b 1 , 

z = (Q-M)z, z(l) = h, 

and, since ce l ^~ N ^b\ = ce*^~ A/ ^6i, it follows that the outputs are the same for all t > 1 as 
well. 

To show that a° and a° are not equivalent, we will show explicitly that the inputs 

{1, 0<t<l, 
0, 1 < s < 1 + s, 
1, t > 1 + 8 

(with varying s > 0) are enough to distinguish the two systems. 

Suppose that the systems have the same output functions under these input functions. 

Let s > and t > 0. Since s(l) = z{\) = h, and x(l + s) = e< Q ~ N \, z(l + s) = e^'^h, 
at time t > 1 + s we have: 

x = Qx + b , x(l + s) = e s( - Q - N h u 
z = Qz + b, z(l + s) = e^-^h . 
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From cx{t) = cz(t) and taking derivatives, one gets 

cQx(t) + cb = cQz(t) + cb \ft>l + s. 

Hence, 

cQx(t) = cQz(t) Vt>l + s. 

Taking more derivatives, we conclude inductively that 

c Q k x(t) = cQ k z(t) Vt > 1 + s, Vk > 0. 

In particular by continuity, we have that cQ k x{\ + s) = cQ k z{\ + s) for all fe > 0, and by 
observability of (Q,c), 

x(l + s) = z{l + s) Vs>0, 

that is, 

e s(Q-N) h = e s(Q-M) bi Vs > Q 

Taking derivatives with respect to s, we conclude that: 

(Q - N) k h = (Q- M)% V k > 0. 

This implies that 

(Q - N)Rx = (Q- M)R 2 , and R x = R 2 , 

where R± = TZ(Q — N,b\), and R 2 = TZ(Q — M,b\). Since (Q — N,b±) is reachable, R\ is 
invertible. From this we conclude that Q — N = Q — M, and hence N = M, a contradiction. 
To complete the proof of Lemma 14.81 we show that a° G T>. First observe that: 

^ = (A + N - M,M,b,c) = (Q-M,M,b,c) 
= MQ,M,b ,c). 

Thus, if we prove that (Q, M,bo,c) £ C, then a° = i/j(Q,M,bo,c) € T>. By Lemma 14.31 
(Q, M, bo, c) € G$. It is thus enough to show that (Q — M, bo, c) is canonical. To see this, note 
that since b\ = det(p(Q))6o and ((Q — N),bo,c) is canonical, it follows that (Q — N,b\,c) is 
canonical. Thus, (Q — M,b\,c) is canonical as (Q — M,bi,c) is similar to {{Q — N)' , d , b^). 
Again, applying the fact that b\ = det(p(Q))bo, one sees that (Q — M,bo,c) is canonical. I 

The above completes the proof of Theorem ^ for the special case r = 1 and a = 1. The 
general case can be obtained by rescaling inputs and time scale, as follows. We consider the 
impulse function u T;a for any fixed r > and a G R. Without loss of generality, we assume 
that For the initial- value problem 

x = (A + aN)x + ab , x(0) = 0, 

let N a = aN, b a = ab, and consider the initial- value problem 

1 



T 



(A + N a )x + b a ), x(0) = 0. 



Then x(t) = x(t/r). It then can be seen that, with 

2V := | (J- A, *N, "b, c) : (A, N, b,c)tV 

and any a £ D a|T , there exists some cr € P QiT such that <r and a have the same output under 
the impulse input u ra , but the two systems are not equivalent. 
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4.3 Proof of Theorem [U 

In this section we consider systems defined as in ©. Let r > and a G R be given. Consider 
the analytic map 

$ : R nxri x R nxn x R nxl x R lxn R nxn x R nxn x R nxl x R lxn 

given by 

(P, N, b ,c)^(P- aN, N, e- rP b , c). 
This is an analytic diffeomorphism whose inverse map is given by 

(A, N, b, c) ^ (A + aN, N, e r{A+aN) b, c). 

Let Ti = &(Qo). The following is a consequence of Lemma 14.21 

The complement Ti c of Ti is the image of an proper algebraic set under a diffeomorphism 
from X to X: Ti c = <&(Gq). Since the image of a measure zero set under a smooth map (see 
e.g. [Q, Lemma 2.6) has measure zero, we conclude: 

Corollary 4.9 The collection Ti is generic. □ 

Lemma 4.10 For any a\ G Ti, there exists 02 G Ti such that 

(a) o\ and 02 have the same output for any u G T Ta \ 

(b) the two systems o\ and 02 are not equivalent. 

Let (Ti = (A, iV, 6, c) G H. Thus, by definition of Ti, there exists (P, N,bo,c) G £0 such 
that o-i = (P - aiV,iV,e- r %,c). Let M = TN'T~ l , where T = T(P,b ,c). Then M ^ N, 
(P, M, 60 , c) G £ , and the system cr 2 := (P - aM, M, e _TP 6 , c) = $(P, M, 6 , c) G ft. 

To prove part (a), pick any u G X TQ! , and assume u(t) = (3 for t > t. The two systems are 
given by 

x = ((P - aN) + uN)x, x(0) = b, 
z = ((P - aM) +uM)z, z(0) = b. 

For < t < t, both systems reduce to the same equation: 

p = Pp, p(0) = b = e~ TP b , 

and in particular, cx(t) = cz(t) = ce tp b for all < t < r. For t > r, the two systems become 

£ = ((P - aiV) + /3iV)x, x(r) = 6 
i = ((P - aM) + /3M)z, z(r) = 6 

and thus have the respective solutions 

z (t) = ce (^+(^-a)M)(t-r) z(r) = ce (P+(/3-a)M)(t-r) bo _ 
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By Lemma 14.31 and the choices of M and N, it follows with 7 = (3 — a that 

which implies that cx(t) = cz(t) for all i > r as well, so part (a) is proved. 

To prove part (b), suppose the two systems a\ and 02 have the same output for all inputs, 
and so, in particular, for all inputs u for which u(t) = a for < t < r. This implies that the 
two systems 

x = ((P — aN) + uN)x, x(t) = e rP b, y = cx, 
z = ((P - aM) + uM)z, z{t) = e rP b, y = cz, 

have the same output for any input. Rewriting the two systems as: 

x = (P + (u — a)N)x, x(t) = e rP b, y = cx, 
z = (P + (u- a)M)z, z(t) = e rP b, y = cz, 

and writing v(t) = u(t) — a for a new input v, one sees that the two systems 

x = (P + vN)x, x{t) =60, y = cx, 
z = (P + vM)z, z(t) =b , y = cz, 

have the same outputs for all inputs v and times > r, which is the same as saying that the 
two systems of type II with associated 4-tuples (P, iV.6o, c) and (P, M, bo, c) are i/o equivalent, 
which is a contradiction in view of Lemma 14.31 and Corollary 14.41 

This completes the proof of Theorem |2 

5 Proofs of Positive Results 

In this section we prove the positive results . 

Let a / be given. Let M. be the set of 4-tuples satisfying the following two properties: 

1. (A, b, c) is canonical, 

2. (A + aN, b) is controllable. 

Since the complement Ai c of M. is defined by the union of the solution sets of the equations 

det n{A, b) = 0, 
det 0(A, c) = 0, 
det K{A + aN, b) = 0, 

we have: 

Lemma 5.1 The complement Ai c of A4 is a proper algebraic set. Consequently, A4 is open, 
dense, and of full measure. 

We will prove that the sets of systems in and S 11 whose 4-tuples are in A4 satisfy the 
conclusions of Theorems 01 and ^respectively. 
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5.1 Proof of Theorem [HI 

Lemma 5.2 Assume that (A, b) is controllable. Then, for any N € M nx?1 and any a 6 



-4, 



is controllable for almost all r > 0. 



Proof. Let A, N, b be given. For any r > 0, 

io £f (fc + 1)! 



where ^(r) = X^fc^o ^it+i)!^ ' ^ Since ^(0) = J, * is analytic, and (A, b) is controllable, 
it follows that (A, \P(t)&) is controllable for almost all r. It then follows that (A, t\P(t)&) is 
controllable for almost all r > 0. I 

To prove Theorem 01 we pick two systems a° = (A±, Nx, b±, c±) and erg = (A 2 , N 2 , b 2 , C2) in 
.A/f, and suppose that they produce have the same output function for each u £ V a . We must 
show that a° = a^. 

Fix a r > 0. Applying u riQ 6 V a to the two systems: 



x 
z 



(Ai + uNi)x + x(0) = 0, y = c\x 

(A 2 + uiY 2 )2 + b 2 u, z(0) =0, y = c 2 z, 



one has: 
where 



cie 



Ai(t-r) 



x(r) = c 2 e 



A 2 (t-r) 



z(r) Vt>r, 



(21) 



(r) = a / e { - M+aN ^ T - s Usb 1 = a e^ Al+aNl >dsh, 



z(t) = a e^ +aN ^ T ~ s Usb 2 = a e^ +aN ^ s dsb 2 . 
Jo Jo 

This holds for any r > 0. 

By Lemma f5.21 we may pick some tq > such that (Ax, x(tq), ci) and (^2, ^(to), c 2 ) are 
both canonical. Since by ()21[) these two triples are i/o equivalent, there exists some invertible 
matrix T € M. nxn such that 



A 2 = T^AxT, z(r ) = T^fa), c 2 = ciT. 
Using in particular that c 2 e A2S = c\e AlS T for all s, Q21|) becomes: 

cie^- r )x(T) = cie^^T^r) Vi > r. 
From the observability of (Ai,c\), it follows that 

x(t) = Tz{t) 



(22) 
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for all r > 0. Equivalently: 

r e{ A 1+a N l)sds bl=T f T e (A 2+ aN 2 )s ds ^ 
JO JO 

for all t > 0. Taking the derivative with respect to r, one gets: 

e (A 1 +aN 1 )r bi = Te {A 2 +aN 2 )T^_ 

Note that this is true for all r > 0. In particular, 

h = Tb 2 . (23) 
On the other hand, taking repeated derivatives in r and then setting r = 0, one obtains: 

(A 1 +aN 1 ) k bi=T(A 2 + aN 2 ) k b 2 VA: > 0. (24) 
This implies, with < k < n — 1, 

U(A 1 + aNi, b x ) = T[U{A 2 + aN 2 , b 2 )], (25) 

and with 1 < k < n, 

(A 1 + aN 1 )[K(A 1 + aN 1 ,b 1 )] = 

T{A 2 + aN 2 )[K{A 2 + aN 2 , b 2 )}. 

Combining this with l)25j) . one sees that 

(A l + aN 1 )[K(A 1 + aNi,b 1 )] = 

T(A 2 + aN 2 )T- l [TZ(A 1 +aN 1 ,b 1 )] 

It then follows from the fact that 1Z(Ai + aN\,bi) is invertible (because {A\ + aN\,b\) is 
controllable) that 

T(A 2 + aN^T- 1 = (A! + aN{). 

It then again follows from (E2J) and the fact that a^0 that N 2 = T^NxT. Combined with 
and (|2Hjl . we have that the systems a° and a 2 are similar, with the similarity matrix given by 
T, and this completes the proof of Theorem |3J 

5.2 Proof of Theorem |U 

The proof of Theorem 0] is almost the same as that of Theorem 01 with Lemma 15 . 21 modified to 
the following: 

Lemma 5.3 Assume that (A, b) is controllable. Then, for any N € W ixn and any a € R, 

'A, e^ A+aN ^b) 

is controllable for almost all r > 0. □ 
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We pick elements o\ = (A±,Ni,bi,ci) and c 2 = (A2, N%, &2j C2) of M. (seen as a class of 
systems of type II) which have the same output function for each u G V a - 

For any r > 0, applying u TjQ , G V a to the two systems: 

x = {A\ + uNi)x, x(0) = b\, y = c\x 
z = (A 2 + uN 2 )z, z(0) = 62, y = C2Z, 

one has, for any r > 0, 

Cl e^ l( *- r) x(r) = c 2 e A2( *- T) z(r) Vi>r, (26) 

where 

x(r) = e^+^Hx, z(r) = e ( A2+aAr2 H 2 . 

By Lemma f5.31 there exists some To > such that both (vli, x(tq), c\) and (A 2 , z(ro), c 2 ) are 
canonical. So, there exists some invertible matrix T G M nxn such that 

= T-^iT, z(r ) = ^^(ro), c 2 = ciT, 

and consequently, (|2*H)) becomes 

cie Al(i - T) x(r) = cie Al( *- T) rz(r) Vt > r. 

For each r > given, using the observability of (Ai, ci), one sees that : 

e (Ai+aATi)T &i _ Te ( J 4 2 +aAf 2 )r &2 

for all r > 0. Starting from here, one can complete the proof by following the same steps as in 
the proof of Theorem |31 

6 A Remark on Sampled Controls 

As remarked earlier, our proofs of the positive results, Theorems 01 and H] were inspired by the 
identification algorithm presented in [T5|. That algorithm aims to find a system equivalent to 
the system being identified, on the basis of observations at discrete instants 0, r, 2r, . . . ,jr, . . ., 
where r is a fixed sampling time, and having applied inputs which have the form Uk T a (f° r 
varying nonnegative integer /c's), i.e., pulses of magnitude a whose width is a multiple of this 
same sampling time r. The motivation is clear: one wishes to use a sample-and-hold strategy, 
which is especially convenient for computer algorithms. Unfortunately, this restriction to fixed 
sampling times means that the algorithm cannot work for generic classes of systems, as we show 
here by means of a counterexample. (Mathematically, the difficulty is that some of the steps 
of the algorithm given in ^21 involve taking logarithms of matrices, which is an ambiguous 
procedure, as the author himself points out on the paper.) To show this shortcoming, for any 
given q > 0, we produce an open class B a of 2-dimensional systems of type I (it is easy to 
generalize to larger dimensions and to systems of type II) with the following properties: for 
every system in B a , there is some other system, which is not i/o equivalent to the original one, 
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yet cannot be distinguished by applying steps of magnitude a and sampled in the above way 
(with fixed r). Thus, our approach, in which r is varied, is actually necessary. 

For any system a° = (A,N,b,c) £ S^, we denote by a° the discrete time system which 
results from sampling the system with r as the length of the sampling interval, and using input 
functions u that are constant over each sampling interval: 

x k+ i = F k x k + u k g k , x(0) = 0, y k = cx k , (27) 

where x k = x{kr),y k = y(kr), u k is the value of u over the interval (kr, (k + l)r) and 

F k = e ^ A+u ^ T , 



9k 



z{ A+u k N)(r~s) hds _ 



For disrete time systems as in ()27(l . the i/o equivalence under a collection of inputs is defined 
in the same manner as in the continuous time case. 

Let B a be the subset consisting of systems of type I for which the 4-tuples 

(A,N,b,c) £l M xl 2x2 xl 2xl xl lx2 

satisfy: 

1. (A + aN, b, c) is canonical; and 

2. A + oiN has a pair of conjugate complex eigenvalues r ± si, with s ^ 0. 

Since the set of 4-tuples for which the triple (A + aN, b, c) is canonical is generic, and the set 
of 4-tuples (^4, N, b, c) for which A + aN has a pair of nonzero conjugate complex eigenvalues 
has a nonempty interior (because of continuity of eigenvalues on matrix entries), the set B a 
contains an open set. 

Proposition 6.1 For any r > 0, and for any a° € B a , there exists some a° € B a that 

1. a° and a° T are i/o equivalent under the collection {u kr ^ a } k >Q, but 

2. a° and a° are not i/o equivalent. 

Proof. Let r > be given, and consider o~° = (A,N,b,c) £ B a . Without loss of generality, 
assume that A + aN is already in "real Jordan canonical form" : 



A + aN 



r —s 
s r 



for some s ^ 0. (If this were not the case, one may simply apply a similarity, and at the end of 
the argument transform back to original coordinates.) Let 



Ao 

Choose an integer / ^ such that 



-2vr/r 
2vr/r 
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• s + 2j2L / 0; and 

• ((A + aN) + ZAo, 6, c) is canonical. 

A generic integer works; note that (since n = 2): 

det 0(A + aN + lA ,c) = det 0(A + aN, c) + det O(lA , c) 

= det 0(A + aN,c) + / det C(A , c). 

and similarly for reachability. Let M be given by 

M = N + -A . 

a 



That is, M is chosen so that 

A + aM = A + aN + /A 



Note that both (A + aN) and (^4 + aM) are is invertible since r + si ^ and r + (s + — )i ^ 0. 
Let £j° = (A,M,b,c), where 6 = (A + aM)(A + aN^b. 

Next we show that the two sampled systems a° and a° T are i/o equivalent under the 
collection {u k T,a}k>o- 

Consider the input function Uk T ,a f° r some k > 0. Clearly the two systems have the same 
state trajectory and same output function if k = (the input is constantly zero). So assume 
k > 1. For j < k — 1, the two sampled systems are given respectively by: 

= Fxj + ag, xo = 0, = ex,, 
Zj+i = FZj + ag, z = 0, yj = cz j: 



where 



F = e ( A + aN ) T = e rT 



cos ST — sm ST 
sin st cos sr 



cos(s + ^)r -sin(s + ^)r 
sin(s + ^)r cos(s + ^?)r 



g(j4+aA/)r = p 



and 



" e (A+«A0(r-9) 6dfl = ( e (A+a7V)r _ J)(^ + a JV)-lft 







( e (A+aM)r_ J ^ i + a¥ |-lj 
e{A+ aM)(T-e)1 de = ? 







(where we have used the fact that e^ T ~ Q ) d6 = (e^ T — I)Q~ l for any invertible Q). Hence, 
Xj = Zj for all < j < k. In particular, outputs coincide at all sampling times t = jr, j < k. 
Over the interval [fcr, oo), the two systems a° and a° are given by 

w = Aw, w(kr) = x(kr) = z(kr), y = cw. 
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This in particular implies that xj = Zj for all j > k and therefore outputs also coincide at 
sampling times t = jr, j > k (as well as for any time t > k). Thus, o° and a° T are i/o 
equivalent under U/% TjCt . 

To show that the two systems a° and a° are not i/o equivalent, consider the constant input 
u = a. Assume the output functions of the two systems are the same under this constant u. 
Then 

c f * e iA+<*N){t-0) had Q = c f' e (A +a M)(t-0) ^ adQ yt > Q) 

Jo io 

and equivalently, 

c (V A+aA ^ - A (A + aN^b = c U A+aM > -A(A + aM^b Vi > 0. 

Since (A + aM) _1 6 = (A + aiV) _1 6, one gets 

c ( e (^+ aAr )* -/)(>! + aiV)- 1 ^ = c (V A + aM )' - (A + aiV)- 1 ^ Vt > 0. 
This in turn implies that 

1 ce (A+«m-e) hd0= j ce (A+aM)(t-9) hdQ y t > 

Jo 

Taking derivatives with respect to t repeatedly, and then setting t = 0, one obtains: 

c{A + aN) j b = c{A + aM) j b Vj>0. 

Since both (A + aN, b, c) and (^4 + aM, b, c) are canonical, it follows that the two systems are 
similar, contradicting the fact that A + aN and A + aM have different pairs of eigenvalues. 

Finally we show that a° £ B a . Note that M was chosen so that (A + aM,c) is observable. 
It is left to show that (A + aM,b) is controllable. Note that P = A + aN and Q = A + aM 
commute. This is because Q = P + — J, and PJ = JP, where 



J 



-1 

1 



Using that PQ = QP and therefore QP" 1 = P _1 Q: 

ft(Q,6) = (6 Q6) = (QP _1 6 QQP^b) = (QP^b QP^Qb) = (QP~ 1 )^(Q,6) 

Since M was chosen so that (Q,b) is controllable, it follows that TZ(Q,b) is non-singular, and 
hence, (A + aM, b) is controllable, as claimed. I 



7 Conclusions and Final Remarks 

For bilinear systems, we showed that step inputs are not enough for identification, nor do single 
pulses suffice, but that that the family of all pulses (of a fixed amplitude but varying widths) 
do suffice. We presented results for single-input single-output systems, since one can obviously 
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identify a multiple-input multiple-output system by considering each pair of input and output 
channels separately, and hence the family of pulses also works for the general case. 

We emphasize that we dealt in this paper with ideal noise-free conditions, and ignored 
stochastic aspects and noisy data, because the underlying theoretical questions of what is ul- 
timately achievable are easiest to understand in a deterministic setting. Tools such as those 
here have been used, however, in the formulation of identification algorithms from noisy data, 
for bilinear systems [22 ■ Nor did we deal here with questions of computational and sample 
complexity. However, the methods used are quite constructive and indeed have appeared in the 
same context in |12j . where numerical implementations are studied; regarding sample complex- 
ity, we leave for further research the generalization of learning-theory results [HI from the 
linear case to the classes of systems considered here. 

Finally, bilinear systems were picked because an elegant result can be established for them, 
as well as their applicability and general interest. However, the study of similar problems to 
those treated here for more general classes of systems is of great interest. 
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